Phase Diagram of Neutron-Proton Condensate in Asymmetric Nuclear Matter 
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We investigate the phase structure of homogeneous and inhomogeneous neutron-proton conden- 
sate in isospin asymmetric nuclear matter. At extremely low nuclear density the condensed matter 
is in homogeneous phase at any temperature, while in general case it is in Larkin-Ovchinnikov-Fulde 
-Ferrell phase at low temperature. In comparison with the homogeneous superfluid, the inhomoge- 
neous superfluid can survive at higher nuclear density and higher isospin asymmetry. 
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It is well-known that the neutron-proton (np) pair- 
ing plays an important role in nuclear physics and as- 
trophysics, such as the structure of medium-mass nu- 
clei produced in radioactive nuclear beam facilities [l[, 
the deuteron formation in intermediate energy heavy-ion 
collisions [3, the pion and kaon condensation the r- 
processjl,!^, and the cooling of neutron stars. The mi- 
croscopic calculations show that the nuclear matter sup- 
ports np Cooper pairing in the "^Si— '^Di channel due to 
the tensor component of the nuclear force, and the pair- 
ing gap is of the order of 10 MeV[i, S i, i, [M [11 [II 
at the saturation nuclear density. At low enough den- 
sity the np Cooper pairs would go over to Bosc-Einstcin 
condensation(BEC) of deuterons in symmetric nuclear 
matter jl,[i]. 

The emergence of isospin asymmetry will generally 
suppress the np pairing, and the condensate will dis- 
appear when the asymmetry becomes sufficiently large. 
Near the saturation density, the np pairing correlation 
depends crucially on the mismatch between the two 
Fermi surfaces, and a small isospin asymmetry can break 
the condensate due to the Pauli blocking effect. At 
very low density, when neutrons and protons start to 
form deuterons and when the spatial separation between 
deuterons and between deuterons and neutrons is large, 
the Pauli blocking loses its efficiency in destroying a np 
condensate. In such situation, the isospin asymmetry can 
be very large, and the np condensate survives in the form 
of deuteron- neutron mixture in momentum space [13, [3] ■ 
Different from the symmetric nuclear matter where the 
thermal motion destroys the np condensate, for asym- 
metric nuclear matter the temperature effect will melt 
the condensate on one hand and increase the overlapping 
between the two effective Fermi surfaces on the other 
hand. As a result of the competition, in a wide density 
regime the temperature dependence of the superfluidity 
is very strange [ill The maximum condensate is not 
located at zero temperature, and the pairing even occurs 
only at intermediate temperature for large isospin asym- 
metry. 

The above results are obtained by assuming the con- 
densate is homogeneous in the ground state. What is the 
true phase structure of np condensate with isospin asym- 
metry when the inhomogeneous Larkin-Ovchinnikov- 



Fulde -Ferrell(LOFF) phasejlSj is taken into account? 
How will the LOFF phase change the strange tempera- 
ture behavior of np condensate found in [1^? In fact, 
there should exist a rich phase structure in asymmet- 
ric nuclear matter, since the isospin asymmetry essen- 
tially plays the same role as the population imbalance in 
the two-component resonantly interacting atomic Fermi 
gas [17]. Different to the cold atoms , in nuclear matter 
the phase separation at large length scale may be forbid- 
den and the LOFF phase may be energetically favored. 
Although the LOFF phase was discussed in asymmet- 
ric nuclear matter at saturation density [T^. the phase 
structure including the LOFF phase in the whole den- 
sity, temperature and isospin asymmetry space and the 
effect of the LOFF phase on the strange temperature be- 
havior of the np condensate are still unknown. We will 
present in this paper the phase diagrams in the density, 
temperature and isospin asymmetry space. 

The often used formulae for superfluid in nuclear mat- 
ter are discussed in detail in [20| where the superfluid 
state is described by a normal and an anomalous nu- 
cleon distribution functions J- and Q. Generally, they 
are functions of momentum p and matrices in spin and 
isospin space. The formulae can be easily generalized to 
study the isospin asymmetric superfluid with total pair 
momentum 2q. We will start with the LOFF phase, and 
the homogeneous phase can be recovered by taking q = 0. 
Like the studies in 

[illil[il,[il[ii, 

we discuss the np 

pairing in the '^Si— '^Di channel with total spin 5* = 1, 
isospin r = and their projections Sz — = 0. In this 
case the distribution functions take the structure 



•^(P) = •?^oo(p)ctoTo + J^o3(p)ctoT3, 

G{P) = ^30(p)a'3CT2T-2, 



(1) 



where ai and are the Pauli matrices in spin and 
isospin spaces. Using the minimum principle of the 
thermodynamic potential and the procedure of block 
diagonalization|2l|, we can express the elements as 

.Foo(p) = l/2-ep[l-./(i?+)-/(i?p)]/(2i?p), 

.Fo3(p) = [/(i?p)-/(i?+)]/2, 

c;3o(p) = -Ap [1 - f{E+) - /(Sp )] / (2i?p) (2) 

with the notations ^p = (p^ + q^) / (2m) — /i, Ep = 



2 



\J Q, + and = Ep ± {Sji + p • q/m), where m is 

the effective nucleon mass in the medium and f{x) = 
l/(e^/"^ + 1) is the Fermi-Dirac function with T be- 
ing the temperature. We have introduced the average 
chemical potential n = {fin + fJ-p)/'^ and the mismatch 
(5/i — {fj,n — /^p)/2 instead of the neutron and proton 
chemical potentials /i„ and /ip. We have also neglected 
the possible neutron-proton mass splitting induced by the 
isospin asymmetry which is believed to be small. The np 
condensate Ap is generally momentum dependent and 
satisfies the gap equation 



A, 



(2^ 



y(p,k)g3o(k), 



(3) 



where V is the nucleon-nucleon (NN) interaction poten- 
tial. The LOFF momentum q should be determined via 
minimizing the the free energy £, which ensures the total 
current js in the ground state to be zero, 



(4) 



where we have used the total nucleon density p = Pn + Pp 
and the isospin density asymmetry Sp — Pn — Pp instead 
of the neutron and proton densities Pn and pp, 



(27r)3 



•?^oo(p), Sp^A 



d^p 
(2^ 



r-^03(p). (5) 



Once the NN potential V is known, we can solve the 
coupled set of gap equations ^ and (g]) together with 
the density equations ([5]) at given temperature T, baryon 
density p or equivalently the Fermi momentum kp = 
(I.Stt^p)^/'^ and isospin asymmetry a — Sp/ p, and obtain 
all possible phases, namely the normal phase Ap — 0, the 
homogeneous superfluid phase Ap 0, q = and the 
LOFF phase Ap 7^ 0, q 7^ 0. By comparing their free 
energies we can determine the true ground state. 

The details of the phase diagram depend on the NN 
potential V we will chose in the numerical calculations, 
however, the qualitative topology structure of the phase 
diagram does not depend on that. To show this, we an- 
alyze the stability of the homogeneous superfluid phase 
against the formation of a nonzero Cooper pair momen- 
tum. For this purpose, we investigate the response of the 
free energy £ to a small pair momentum q via the small 
q expansion, £{q) — S{0) + js ■ c[/m + psq^/{2m) 
where — md£ /d(\ is the total current which is pro- 
portional to the left hand side of (|4]), and ps is just the 
superfluid density defined by ps = md^E/dc^ of which 
the explicit form reads 



d3 



P P 



m J (27r)3 3 



[/'(i?p+) + /'(i?p)]|_ (6) 



with the definition f'{x) = df{x)/dx. The current 
vanishes due to the gap equation for q, and the sign of 
Ps controls the stability of the homogeneous superfluid. 



i.e., a negative ps means the LOFF phase has lower free 
energy than the homogeneous superfluid phase. 

The momentum dependence of the gap function Ap is 
normally rather weak in a wide momentum region[l3|, 
and we can approximately treat it as a constant A for 
a qualitative analysis. At zero temperature, once the 
isospin asymmetry is turned on, there must exist a sharp 
breached region where the quasiparticle energy E~ < 
with the necessary condition Sp > A, and the momentum 
integration in ^ can be analytically carried out, 



Ps= P 



1 - 



Sp 



Stt^P 



, (7) 



where p± = ^^2mp± 
p± = p± ^5p^ - A2 



are possible gapless nodes with 
At high density the matter is in 



BCS regime where Sp, A ^ p, and the breached region is 
P- < IpI < p+. Since p± are close to the Fermi momen- 
tum kp, the superfluid density should be negative since 
Ps ~ p{l — Sp/ y/Sp"^ — A^). On the other hand, at low 
enough density the matter is in BEG regime, the chemical 
potential p becomes negative which leads to p- < and a 
reduced breached region < |p| < p+- In this case p+ is 
much smaller than kp and the superfluid density becomes 
positive. Therefore, at zero temperature the superfluid is 
expected to evolve from an inhomogeneous phase to the 
homogeneous phase when the nuclear density decreases, 
which is a general phenomenon for BCS-BEC crossover 
with population imbalance [l7j. 

At finite temperature, the breached region is smeared 
due to the thermal excitation, and Sp > A is not nec- 
essary. At the critical temperature Tc there should be 
a second order phase transition from the superfluid to 
normal state, and for temperature T < Tc the pairing 
gap behaves as A(r) cx (1 - T/TcY^'^ which leads to 
the regular behavior of the superfluid density Ps{T) cx 
(1 — T/Tc) > 0. This means the temperature tends to 
stabilize the homogeneous phase. Combining with the 
behavior of the superfluid density at zero temperature, 
the homogeneous phase at sufhciently low density will 
keep stable at any temperature below Tc, while at high 
density there must exist a turning temperature Ts where 
Ps changes sign, and the superfluid should be in inhomo- 
geneous phase at low temperature T < Ts and in homo- 
geneous phase at high temperature < T < Tc- Since 
the homogeneous state is unstable at low temperature, 
combining with the fact that the critical isospin asym- 
metry for the LOFF phase is much larger than the ho- 
mogeneous phase ^195, the strange temperature behavior 
of the pairing gap found in ^15i] is probably unrealistic. 

We now move to numerical calculations. The Paris NN 
potential is often used to describe the nuclear structure 
and nucleon superfluidity, and describes well the BCS- 
BEC crossover of np condensate^, [l^. Since the quali- 
tative topology structure of the phase diagram does not 
depend on specific models, for the sake of simplicity, we 
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FIG. 1: The pairing gap A for homogeneous (dashed lines) 
and inhomogeneous (sohd hnes) condensates and LOFF mo- 
mentum q as functions of temperature at normal density po 
and for two values of isospin asymmetry. 



employ a separable form of the Paris NN potential 

7n 



y(ri,r2) = Wo 



1—77 




5{vi-V2), (8) 



which was developed in (22f] to reproduce the pairing gap 
in 5 = l,r = channel and the bound state between 
zero energy and deuteron binding energy, where po is the 
normal nuclear density and the parameters VQ,r],j and 
an energy cutoff Sc to regulate the model are determined 
by recovering the pairing gap in the realistic Paris NN 
potential. This separable form can also describe well 
the BCS-BEC crossover of np condensate 14]. In our 
numerical calculation, we choose vq = —530 MeV-fm^, 
77 = 0, ec = 60 MeV and take m as the density-dependent 
nucleon mass corresponding to the Gogny force DlSf23|. 
We have checked that different parameter setsi22:i lead to 
only a slight change in the numerical results. 

In the low density BEG regime, the homogeneous phase 
is stable at any temperature below Tc, and the conden- 
sate is a regular decreasing function of temperature. Be- 
yond this regime the strange phenomenon of the homo- 
geneous condensate arises: The maximum pairing gap 
is located at non-zero temperature, and the condensate 
even appears only at intermediate temperature with two 
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FIG. 2: The phase diagrams in kp — a, kp — T and a — T 
planes. In each plane, the labels HS, LOFF and Normal in- 
dicate homogeneous superfluid, LOFF superfluid and normal 
phase, the dashed line is the border of the unstable HS, and 
the three phases meet at a Lifshitz point L. 



critical temperatures To and Tc for large isospin asym- 
metry, which is shown as dashed lines in the upper panel 
of Fig[l]for saturation density without loss of generality. 
However, the superfluid density ps for the homogeneous 
phase is negative at T < Tg and positive at < T < 
where the turning temperature Ts is larger than the lower 
critical temperature Tq. This tells us that the homoge- 
neous phase is stable at high temperature T > Tg and 
unstable to formation of LOFF condensate at low tem- 
perature T < Ts- By calculating the LOFF pairing gap 
A and momentum (? = |q| which are shown as solid lines 
in the upper and lower panels of FigU] and comparing 
the free energies for the homogeneous and LOFF states, 
the LOFF phase is energetically more favored than the 
homogeneous phase a.t T < Tg. Especially, different to 
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the homogeneous condensate, the LOFF condensate al- 
ways starts at zero temperature. Therefore, after con- 
sidering both homogeneous and inhomogeneous conden- 
sates, the strange temperature dependence of the pairing 
gap[l^[l^ disappears and the condensate becomes a reg- 
ular decreasing function of temperature. The LOFF mo- 
mentum q, shown in the lower panel of FiglU drops down 
with increasing temperature and approaches to zero con- 
tinuously at Ts , which indicates a continuous phase tran- 
sition from homogeneous phase to LOFF phase. The 
continuity can be proven analytically. The gap equation 
^ for q can be written as qW{q) = with a trivial so- 
lution q — for the homogeneous phase and a non-zero 
solution from W{q) — for the LOFF phase. Using the 
expansion for £ we find ps — 14^(0). Therefore, we must 
have q = a,t T = Ts, providing that the LOFF solution 
is unique. 

The phase diagrams of the np pairing are shown in 
Figl21 We first discuss the one in the kp — a plane at a 
very low temperature T = 0.1 McV. When p ^ we find 
/i^— efc/2ata = where ej, is the deuteron binding 
energy, which means that the np condensate survives in 
the form of deuteron BEG. Consistent with the findings 
in atomic Fermi gas|17|, the homogeneous phase (HS) 
is stable only at very low density BEG regime. In this 
regime, the critical isospin asymmetry can be very large, 
and even approaches to 1 for kp < 0.23 fm~^. Beyond 
this extremely low density regime, the superfluid density 
of the homogeneous phase becomes negative which indi- 
cates that the LOFF phase is energetically favored. By 
calculating the superfluid density of the HS phase and 
the LOFF solution, we can determine the phase bound- 
aries between HS and LOFF and between LOFF and 
normal phase. The LOFF momentum is large at high 
p and high a and approaches to zero at the HS-LOFF 
boundary which means a continuous phase transition. If 



we consider HS only, the HS-Normal boundary (dashed 
line) is below the LOFF-Normal boundary, this shows 
that the introduction of LOFF phase enlarges the super- 
fluid region. In the kp — T and a — T planes, the HS and 
LOFF phases are separated by the turning temperature 
Tg. Note that starts at kp ^ in kp ~ T plane which 
corresponds to the stable HS at extremely low density 
but starts at a = in a — T plane for high density which 
means that the HS with small isospin asymmetry is easy 
to be stabilized. Again, in comparison with the calcula- 
tion with only HS, the superfluid is extended to higher 
density or higher asymmetry region due to the introduc- 
tion of the LOFF phase, and the unstable HS-Normal 
boundary (dashed lines) which reflects the strange "in- 
termediate temperature superfluidity" is replaced by the 
LOFF-Normal boundary. The phase transitions in the 
three planes are all of second order, and in any case the 
three phases meet at a Lifshitz point L 24]. The a — T 
phase diagram we obtained is very similar to the generic 
phase diagram of two-component ultracold Fermi gas in 
a potential trap [2^. 

In summary, we have qualitatively investigated the 
phase structure of np condensate in isospin asymmetric 
nuclear matter and confirmed our analysis with a model 
NN potential. The important findings are: 1) The LOFF 
phase is the ground state in a wide region of nuclear den- 
sity, temperature and isospin asymmetry, except for very 
low density and high temperature. 2) The strange tem- 
perature behavior of the np condensate is washed out by 
the LOFF phase at low temperatures. 3)The superfluid 
region is expanded to high density and high asymmetry 
due to the introduction of LOFF phase. 
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